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statement  of  Problem 


The  knowledge  of  the  dynamics  of  elementary  chemical  reactions  occurring  in  the 
gas  phase  and  on  surfaces  is  crucial  to  the  understanding  of  the  more  complex  mechanisms 
of  decomposition  of  energetic  materials.  Heterogeneous  reactions  control  the  initial 
decomposition  of  solid  propellants  and  explosives  and  can  also  play  a  role  in  combustion  in 
closed  environments  in  which  wall-sensitive  reactions  are  important.  Gas-phase  reactions 
in  the  primary  flame  zone  immediately  above  the  surface  of  the  propellants  and  explosives 
control  the  heat  conduction  back  to  the  solid  material  and  thereby  influence  the  stability  and 
efficiency  of  the  energetic  materials.  Exothermic  gas-phase  reactions  of  initial 
decomposition  products  release  heat  which  is  then  absorbed  by  the  solid,  thus  determining 
the  buming-rate-versus-pressure  characteristics  of  the  combustion. 

Major  limitations  in  the  understanding  of  the  complex  combustion  mechanisms  of 
the  decomposition  of  energetic  materials  are  the  lack  of  knowledge  of  all  the  important 
reaction  pathways,  the  lack  of  knowledge  of  dynamical  constraints  of  individual  elementary 
reactions,  and  the  lack  of  accurate  kinetic  data  of  the  elementary  reaction  steps.  The  goal  of 
the  research  reported  here  is  the  development  of  theoretical  methods  to  aid  in  the 
interpretation  and  modeling  of  gas-phase  combustion  processes.  The  methods  are  based 
upon  variational  transition  state  theory  (VTST),  which  has  been  extensively  tested  for  gas- 
phase  reaction  involving  only  a  few  atoms.  The  main  thrusts  of  this  research  have  been 
continued  tests  and  development  of  the  theory  with  the  intent  of  extending  the  regions  of 
applicability  of  the  methods,  extensions  of  the  methods  to  treat  reactions  involving  larger 
polyatomic  species,  and  applications  to  systems  of  practical  importance.  The  study  of 
heterogeneous  reactions  has  required  further  extensions  of  the  methods.  The  aim  in  the 
study  of  heterogeneous  reactions  has  been  the  qualitative  understanding  of  the  dynamics  of 
reactions  occurring  on  well  characterized  crystallographic  surfaces  and  the  systematic  study 
of  qualitative  features  of  these  processes. 


Summary  of  Important  Results 


Exciting  new  advances  have  been  made  within  this  three  year  research  program 
which  have  brought  us  closer  to  our  goal  of  a  practical,  yet  accurate  method  for  computing 
rate  constants  for  a  wide  variety  of  chemical  reactions.  This  work  is  roughly  divided  into 
continuing  tests  of  the  accuracy  of  the  methods  of  variational  transition  state  theory 
(VTST),  new  theoretical  developments  and  numerical  methods,  and  applications  to  systems 
of  fundamental  importance. 

Quantitative  tests  of  VTST. 

Although  many  quantitative  comparisons  have  previously  been  made  between  rate 
constants  computed  using  VTST  and  those  obtained  from  accurate  quantum  mechanical 
calculations,  most  of  these  have  been  for  model  collinear  reactions.  With  the  advent  of  the 
new  generation  of  supercomputers  and  the  development  of  new  methods  of  quantum 
scattering  theory,  there  has  been  a  recent  flourish  of  activity  in  the  area  of  accurate  quantum 
calculations  for  realistic  reactions  in  three  dimensions.  New  quantum  mechanical  results 
for  several  isotopic  variants  of  the  reactions  H  +  H2  and  O  +  H2  have  allowed  further  tests 

of  our  VTST  calculations.  These  tests  have  once  again  confirmed  the  utility  of  VTST. 
Both  sets  of  reactions  have  high  classical  barriers  and  are  dominated  by  quantum 
mechanical  tunneling  at  low  temperatures.  Even  so,  the  methods  were  shown  to  be 
accurate  to  within  about  40%  even  down  to  room  temperature. 

Perhaps  the  most  interesting  tests  of  VTST  have  been  comparisons  with  reliable 
experimental  rate  constants.  Such  comparisons  require  accurate  potential  energy  surfaces 
for  carrying  out  the  theoretical  calculations.  For  the  reaction  of  0(^P)  with  H2  such  a 

potential  exists  and  VTST  calculations  have  been  carried  out  for  the  rate  constants  and 
deuterium  kinetic  isotope  effects  over  a  temperature  range  extending  from  below  room 
temperature  up  to  2400K.  Since  then  definitive  experimental  studies  have  been  carried  for 
this  reaction  measuring  the  thermal  rate  constant  at  low  temperatures  and  also  directly 
measuring  the  deuterium  kinetic  isotope  effects.  These  measurements  combined  with 
recent  shock  tube  measurements  at  higher  temperatures  have  given  a  complete 
determination  of  the  rates  over  an  extremely  wide  temperature  range.  The  comparison  of 
the  measured  and  theoretical  rate  constants  have  indicated  the  predictive  powers  of  the 
theory.  However,  an  even  more  exciting  conclusion  from  these  studies  is  the  experimental 
confirmation  that  quantum  mechanical  tunneling  plays  a  dominate  role  in  this  reaction  at 
temperatures  from  room  temperature  to  about  600K. 
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Hydrogen  transfer  reactions  are  typically  facile  and  are  ubiquitous  in  combustion  of 
hydrocarbons  and  most  conventional  propellants  and  explosives.  For  the  transfer  of  H 
between  two  heavy  moieties,  the  angle  between  the  bottom  of  the  reactant  valley  and  the 
bottom  of  the  product  valley  in  a  mass  weighted  coordinate  systems  is  small.  A 
consequence  of  this  small  skew  angle  is  that  the  reaction  path  must  go  through  a  highly 
curved  region  in  proceeding  from  reactants  to  products.  We  have  found  that  the  effect  of 
variationally  optimizing  the  location  of  the  geno-alized  transition  state  dividing  surface  (the 
dynamical  bottleneck  for  the  reaction)  can  be  very  important  for  such  small  skew-angle, 
high  reaction-path-curvature  systems.  The  importance  of  the  variational  procedure  has  not 
been  definitively  tested  in  previous  studies  because  in  those  studies  the  systems  have  had 
classical  barriers  which  were  high  enough  that  quantum  mechanical  tunneling  dominated 
the  reaction  over  a  wide  temperature  range.  The  reaction  I  +  HI  has  been  studied  on  a 
potential  energy  surface  with  a  low  classical  barrier.  For  this  reaction  tunneling  effects 
were  found  to  be  nominal  and  comparison  with  accurate  quantum  mechanical  results  have 
shown  that  the  dynamical  bottleneck  must  be  positioned  away  from  the  symmetrically 
located  saddle  point  to  obtain  quantitative  agreement  between  the  quantum  mechanical  and 
VTST  results. 

New  theoretical  and  numerical  developments 

For  reactions  with  moderate-to-large  reaction  path  curvature  in  which  quantum 
mechanical  tunneling  is  important,  the  least-action  ground-state  (LAG)  method  was 
previously  shown  to  provide  insight  into  the  dynamics  of  the  reaction  as  well  as  an  accurate 
method  for  computing  a  quantum  mechanical  tunneling  correction  factor.  This  method  was 
previously  tested  for  reactions  which  were  nearly  thermoneutral  and  for  which  the 

V 

correction  factor  is  dominated  by  contributions  from  tunneling  from  the  ground  vibrational 
state  of  reactants  into  the  ground  vibrational  state  of  products.  For  exothermic  reactions 
such  as  Cl  +  HBr  there  is  a  propensity  to  populate  excited  vibrational  states  of  the  products. 
The  LA  method  was  extended  to  treat  this  inherently  vibrationally  nonadiabatic  process  of 
tunneling  into  vibrationally  excited  states  and  successfully  tested  by  comparison  with 
accurate  quantum  mechanical  calculations  for  the  model  collinear  reaction  Cl  +  HBr  -¥  HCl 
+  Br. 

Phenomenological  models  based  upon  conventional  transition  state  theory  with 
approximate  tunneling  correction  factor  have  been  widely  used  over  the  last  40  years  to 
interpret  kinetic  isotope  effects  in  hydride  transfer  reactions.  The  transfer  of  a  light  species 
between  two  heavy  moieties  manifests  small  skew  angles  and  large  reaction  path  curvature. 
The  insight  into  such  processes  gained  from  our  VTST  studies  including  the  LAG 
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tunneling  corrections  have  shown  that  these  phenomenological  models  are  based  upon 
incorrect  dynamical  assumptions.  However,  we  have  shown  how  variational  effects  and 
large-curvature  tunneling  effects  can  be  incorporated  into  a  phenomenological  model  of 
kinetic  isotope  effects  to  give  a  more  consistent  picture  of  the  underlying  dynamics. 

The  calculadon  of  thermal  rate  constants  samples  collision  energies  which  are  near 
to  the  maximum  of  the  ground-state  adiabatic  barrier.  Reactions  involving  vibrationally 
excited  reactants  sample  higher  collision  energies  and  provide  a  probe  of  wider  regions  of 
the  potential  energy  surface.  In  addition,  these  reactions  provide  a  more  stringent  tests  of 
the  VTST  methods.  Approximate  methods  for  treating  excited  state  reactions  have  been 
developed  based  upon  VTST  and  tested  for  several  model  and  realistic  reactions.  It  was 
found  that  an  adiabatic  description  of  the  reaction  in  the  reactant  channel  up  to  the  first 
occurrence  of  large  reaction-path  curvature  provided  a  convenient  picture  of  the  dynamics 
which  is  capable  of  semiquantitative  accuracy. 

Another  area  of  crucial  interest  is  the  application  of  VTST  to  reactions  involving 
larger  polyatomic  species  which  will  be  important  in  combustion  processes.  For  these 
larger  reaction  systems  global  potential  energy  surfaces  will  generally  not  be  available  and 
methods  were  previously  developed  to  allow  accurate  ab  initio  electronic  structure 
information  to  be  directly  used  in  the  VTST  calculations.  Methods  for  interpolating  the  ab 
initio  information  within  the  framework  of  the  reaction  path  Hamiltonian  (RPH)  on  a 
sparse  grid  allows  the  calculation  of  accurate  rates  from  only  a  nxxlerate  amount  of  ab  initio 
information.  Quantitative  tests  of  these  RPH  interpolation  methods  have  only  recently  been 
carried  out  by  comparing  rate  constants  computed  using  a  global  potential  energy  surface 
with  those  computed  using  RPH  interpolation  of  pseudo-a&  initio  information  obtained 
from  the  global  potential  energy  surface.  These  calculations  indicate  that  rate  constants  for 
the  four  atom  reaction  OH  +  H2  can  be  obtained  to  a  numerical  accuracy  of  about  25%  at 

temperatures  down  to  room  temperature  using  information  about  the  potential  and  its  first 

t 

and  second  derivatives  at  8-10  points  along  the  reaction  coordinate. 

Another  concern  in  the  use  of  RPH  interpolation  techniques  is  how  accurate  the 
reaction  path  (the  minimum  energy  path  in  a  mass-weighted  coordinate  system)  is  located. 
This  is  determined  by  the  numerical  integration  scheme  used  to  follow  the  path  of  steepest 
descents  and  the  step  size.  We  have  made  a  comparison  of  several  numerical  integration 
algorithms  and  have  found  that  to  obtain  a  numerical  accuracy  of  15%  the  simplest  fixed- 
step  method  allowed  the  largest  steps  to  be  taken.  However,  for  much  higher  accuracy, 
0.01-0.1%,  a  variable-step  predictor-corrector  algorithm  was  shown  to  be  the  most 
efficient.  These  algorithm  were  tested  for  the  reactions  OH  +  H2  and  CH3  +  H2  using 
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global  potential  energy  surfaces  and  for  the  reactions  CH3  +  Hj  and  HCN  ->  HNC  using 
ab  initio  electronic  structure  information 

The  VTST  methods  have  been  extended  to  allow  treatment  of  reactions  occurring  on 
surfaces.  This  includes  simple  processes  such  as  hydrogen  migration  on  surfaces  and  the 
more  dynamically  complex  process  of  dissociative  chemisorption.  A  major  advantage  of 
VTST  over  other  dynamical  methods  for  treating  these  systems  is  that  the  motion  of  the 
surface  atoms  can  be  easily  included  in  the  calculation,  thereby  incorporating  the  effect  of 
phonons  on  the  reaction  process. 

Applicatigns 

An  exciting  prospect  for  obtaining  reliable  rate  constants  for  systems  of  practical 
importance  in  the  study  of  energetic  materials  has  been  collaboration  with  quantum 
chemists  to  obtain  ab  initio  electronic  structure  information  as  input  into  the  VTST 
calculations  Such  a  program  is  underway;  in  collaboration  with  Drs.  Byron  Lengsfield 
and  George  Adams  at  the  Ballistic  Research  Laboratory  at  Aberdeen  Provings  Grounds  we 
are  currently  studying  the  isomerization  reaction  HONO  ->  HNO2.  HONO  is  known  to  be 

an  important  species  in  the  decomposition  of  propellants  and  explosives  containing  nitro 
substituents.  Cunently  the  mechanism  of  depletion  of  HONO  is  not  well  understood  and  it 
has  been  suggested  that  the  isomerization  to  form  HNO2  may  be  feasible.  Preliminary 

indications  from  our  studies  indicate  that  the  reaction  rate  is  extremely  small  below  about 
600K  but  may  be  nonnegligible  at  higher  temperatures.  Because  of  the  high  classical 
barrier,  even  at  600K  the  reaction  is  dominated  by  quantum  mechanical  tunneling. 

Progress  has  also  been  made  in  the  study  of  reactions  occurring  on  surfaces.  A 
major  difficulty  in  these  studies  is  the  lack  of  accurate  information  about  the  potential 
energy  surface.  Because  of  the  extremely  large  number  of  atoms  in  the  system  the 
prospects  of  obtaining  accurate  ab  initio  electronic  structure  information  in  the  near  future 
is  bleak.  The  most  promising  approach  currently  is  provided  by  semiempirical  potentials 
which  are  designed  to  have  the  correct  topology  and  are  empirically  fit.  This  is  similar  in 
spirit  to  the  approach  taken  to  gas-phase  reactions  twenty  years  ago  in  which  LEPS 
potentials  were  empirically  adjusted  to  reproduce  experimental  results.  These  potential 
surfaces  could  then  be  used  in  theoretical  calculations  of  other  dynamical  quantities  and  to 
provided  insight  into  the  dynamics  of  the  reaction.  A  particularly  fruitful  approach  to 
treating  the  bulk  properties  of  metal  surfaces  has  been  the  use  of  the  embedded  atom 
method  (EAM)  which  is  based  upon  density  functional  theory.  This  method  has  also  been 
extended  by  Baskes,  Daw,  and  Foiles  to  treat  the  interaction  of  adsorbates  such  as 
hydrogen  with  a  metal  surface.  As  a  test  of  the  accuracy  of  the  .semiempirical  potential  for 


H  interacting  with  bulk  Ni,  we  have  studied  diffusion  of  H  on  the  100  and  1 1 1  surfaces  of 
Ni  and  diffusion  in  bulk  Ni.  The  computed  rates  are  in  quantitative  agreement  with  the 
plentiful  experimental  measurements  of  bulk  diffusion.  For  diffusion  on  the  100  surface 
the  experimental  measurements  are  much  harder  and  less  reliable,  but  our  calculated  results 
are  in  qualitative  agreement. 

The  description  of  dissociative  chemisorption  such  as  H2  on  Ni(lOO)  is  complicated 
by  the  need  to  describe  the  covalent  H2  bond  in  the  gas  phase  and  the  H  interactions  with 
Ni  when  it  is  chemisorbed.  Covalent  bonding  is  not  easily  described  within  density 
functional  theory  and  therefore  the  EAM  method  by  itself  is  not  adequate  for  describing  the 
potential.  We  have  developed  a  general  method  of  incorporating  the  EAM  description  of 
the  metal-metal  and  metal-adsorbate  interactions  within  a  diatomics-in-molecules  (DIM) 
framework  that  allows  a  flexible  framework  for  semiempirically  fitting  the  gas-surface 
potential.  We  note  that  the  DIM  model  is  equivalent  to  the  LEPS  potential  when  only  two 
diatomic  potentials  are  used  to  describe  the  atom-atom  interactions.  Thus,  the  embedded 
DIM  method  is  just  a  natural  extension  of  LEPS  to  treat  gas-surface  interactions.  Work  is 
now  in  progress  to  study  H2  chemisorbing  on  Ni(lOO) 
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